A Cloudy Bag Model for the 
S-D Wave K- N System 



Guangliang He and Rubin. H. Landau 
Physics Department, Oregon State University, Corvallis, OR 97331 
(February 9, 2008) 

The cloudy quark bag model is applied to the coupled {KN, Ett, Att) system for S — D partial 
waves. Energy-dependent separable potentials are derived, as are new numerical algorithms for 
solving the coupled Lippmann-Schwinger equations. The parameters of the model are fit to K^p 
scattering and reaction cross sections, branching ratios, and mass spectra from K~p Etttttt, Atttttt. 
Within the constraints of the model, the branching ratio and Ett mass spectrum data are in conflict. 
The Epi3(1385) and A_do3(1520) resonances are found to be predominately elementary bag states 
with considerable dressing for the Epia. The Asoi(1405) appears as a complicated composite systems 
• arising from two poles. The model with certain parameter sets does predict two sign changes in 

0^ ' the real part of the KN scattering amplitude near threshold, but they are not quite at the correct 

0^ , energies to produce agreement with the sign of the strong interaction shift of kaonic hydrogen. 

^ . I. INTRODUCTION 

<: 

It is generally accepted that hadrons are composed of quarks and that strong interactions among quarks will 
[ be explained by quantum chromodynamics (QCD). Unfortunately, QCD is still too complicated to solve and so 
' we construct phenomcnological models which incorporate confinement, asymptotic freedom, and the symmetries of 
J> QCD — yet still permit analysis of experimental data and application to nuclear few-body systems. The MIT "bag" 
is a model for a bubble in the QCD vacuum and imposes confinement by placing the quarks in a finite bag of 
about 1 fm radius. If the model were interpreted literally, this large a bag would shatter conventional nuclear physics 
since nuclei would be filled more with quarks than nucleons. 
' Soon after its introduction it was realized that the bag model violates chiral symmetry (conservation of quark 
[ helicity) — the second best symmetry (^ 7%) of the strong interactions ||^. The violation arises from quark reflections 
^\ off the bag wall changing the quark's momentum but not its spin. The symmetry is restored by coupling in a pion field 
_ to the bag's surface In the "little brown bag" of Brown and Rho |Q, an interior phase contains free, massless 

. quarks, while an exterior phase contains pions but no quarks. The pressure of the pion field on the outside of the 
,_L| ' bag then compresses the bag down to the very small 0.3 fm little brown bag. Clearly with this small a bag there will 
O : be very few quark effects expected inside of nuclei, and the conventional views of nuclei as collections of nucleons 
exchanging mesons holds. 

In the cloudy bag model (CBM) of Miller, Thomas, and collaborators [^,^, the exterior pion field gets quantized 
and penetrates the bag. While the CBM's Lagrangian is non-linear, practical calculations are carried out as a 
perturbation expansion with the zeroth order term yielding the MIT bag model. The first tests of the CBM dealt 
with static properties of hadrons, such as charge radii and magnetic moments, and were successful 0,^. Latter 
' tests dealt with dynamic properties, such as pion-nucleon scattering, and required an extension of the Lagrangian 
to incorporate volume in addition to surface coupling M. Since the extension incorporated Weinberg's effective ttN 
Lagrangian, it is not surprising that success was found 

In order to study the KN and KN systems, strangeness was incorporated into the CBM by extending the Lagrangian 
from SU(2) to SU(3) |12 1^. However, the complications of thresholds and exotic KN resonances led Veit et al. 
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to conclude that the model "has trouble with resonances" for both the KN and KN systems. For example, they 
found a KN scattering amplitude which appeared to have a A resonance signal approximately 5 MeV below threshold 
rather than the expected 27 MeV below. Further analysis by Fink et al. indicated that the resonance signal arises 
from the influence of the threshold cusp on a pole some 13 MeV above threshold. 

The KN findings are interesting because they reflect on the nature of the subthreshold A resonance and possibly 
on the puzzling experimental measurements of the IS* strong interaction level shift in kaonic hydrogen p^-p^]. The 
hydrogen experiments indicate that the real part of the K^p scattering amplitude has either no or two sign changes 
right above threshold ]T5| , p^ -pl| , whereas K matrix analyses and most potential models predict one sign change. The 
CBM fits have two sign changes close to threshold, but on either side of threshold. 

It has not been clear to us what level of agreement to expect between models of the KN coupled system and 
energy-dependent scattering and reaction data. On the one hand the CBM assumption of a rigid spherical bag and 
of point-like pion fields is possibly too simple, yet on the other hand it may be that only more partial waves are 
needed to get the energy dependence right. There is, after all, experimental indication that significant P wave enters 
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at 184 MeV/c and there is the D-wave resonance at 390 MeV/c. To describe additional data and to test the 
model further, we have extended the KN cloudy bag model to include the P and D partial waves. This extends the 
contact interaction to include space-like terms in addition to time-like ones, and the general interaction to include 
spin-orbit forces. We have deduced effective potentials for use in a Lippmann-Schwinger equation, and have fit the 
potentials' parameters to data up to 520 MeV/c. We were aided in our work by the similar kaon work of Veit et al. 
||13| ; similar, yet with the additional and stronger interactions available to the antikaon making for more differences 
than similarities. 



B = --n-d[q{x)q{x)\^^^ (4) 



II. THE CLOUDY BAG LAGRANGIAN 

The MIT bag model's Lagrangian is 

i _^ 1 _ 

Cmit{x) = i-q ^ q - B)dy ~ -qqSs (1) 

Here q is the quark field, B is a universal constant called the "bag pressure", and for a nondeformable, static, 
spherical bag, the bag ansatz is imposed by setting the volume step-function 9y = 9{R — r) and the surface delta 
function Sg — 6{R — r). By demanding that the action S = Jd'^x C{x) remains stationary under arbitrary changes in 
5, q, 9y, and (5^, we obtain the equations of motion: 

I ^q{x) = 0, (r < i?), (2) 
i ■/iq{x) = q{x), {r ^ R), (3) 
1 

—1 
2 

For the static case, the unit normal to the bag's surface is n™" — (0, f). Equation ^ is the Dirac equation inside the 
bag, is the linear boundary condition of confinement, and is the stability condition on the bag pressure. We 
shall use solutions to these equations in constructing the baryon-meson effective potentials. 

The Lagrangian (|l|) does not preserve chiral symmetry since the surface term ^qqSg reverses the helicity of the quark 

hitting the bag's surface. To incorporate chiral symmetry, a pion field (p is coupled throughout the bag's volume [^^l): 

CcBM ^(^y^q~B^ 0y - lqe^^-$-f'^/fq6, + 1 (^D^$) ' (5) 

Here the arrow over (j> indicates its isovector nature, r is a vector of Pauli matrices for SU (2) (isospin), / is the meson 
octet decay constant, and is the covariant derivative: 

= id^^)4> + f sinW f)d^4> (6) 

1/2 

■ ^/^ (7) 



{$■$)' 



While the Lagrangian has no obvious predictions for low-energy pion-nucleon scattering, Thomas transformed 
it to one containing a volume coupling of the covariant derivative of the quark field: 

CcBM = C-qpq~ B)0y - ^qqS, + ^ (i?^<^) ' + ^f^^'^^^q ■ {0^)0, (8) 



D^,q = d^q 



cos(0//) - 1 



(0 X 9^0) q (9) 



2 

This Lagrangian incorporates the Weinberg- Tomozawa relation for zero-energy S'-wave pion scattering. 



III. APPLICATION TO KN SCATTERING 

To apply the Lagrangian (^) to KN scattering, Veit et al. [|2[jl4j extended the internal symmetry from flavor 
SU{2) X SU{2) to SU{3) x SU{3). The quark field q then becomes any member of the SU(3) triplet (u,d,s), the 
meson field (j> any member of the octet, the SU{2) Pauli matrices f are replaced by SU{3) Gell-Mann matrices A, and 
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the cross product in the covariant derivative includes SU{3) structure constants {A x B)c = J2a b fabcA-aBi,. Since 

(f = yields the MIT bag Lagrangian (|^), which was successful for static baryon properties, Veit ct al. |^ postulated 
that if the energy is low then small (f> should be a good approximation for scattering. Accordingly, they expanded 
^CBM around (p — to obtain the linearized, volume-coupled, SU(3)xSU(3), CBM Lagrangian: 
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CBM 



= c 



MIT 
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^M = ^{d,,(, 



(10) 
(11) 
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where Lmit (|l|) describes the free bag, Cm the free meson field, Cg the s-channel interaction, and Cc the contact 
interaction. 

In our application of the CBM to the KN system, we restrict the energy to (1250 < E^m < 1550)MeV, that is, to 
both sides of the KN threshold at 1432 MeV. For these energies we consider the six, strangeness —1, baryon-meson 
{BM) channels which couple strongly: 



K-p 



K-p {threshold") 1432 MeV 

K^n -5 MeV 

E-7r+ +95 MeV 

E^ttO +104 MeV 

E+TT- +103 MeV 

Att*' +181 MeV 
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This system supports a number of resonances which appear in different channels below and above threshold. For low 
energies the relevant ones are Al/2j=soi(1405) and Epi3(1385) below threshold and Ado3(1520) above. 

The calculations to follow get rather complicated and so we have placed many of the details in appendices. We will 
describe how we extended the CBM up to D waves and in each partial wave deduced effective potentials to be used 
in the Lippmann-Schwinger equation 



r^4k',/i';k,Ai) =F^a(k',A*';k,^)+;^ d^p 



F^-y (k', ii'; p, y)T^a{p, v\ k, p) 
E + ie- E^ (p) 



(15) 



Here the incident and final particles are in channels a and (3 (defined in table ||), have spin and ^' , and COM 
momentum k and k'. The intermediate state is channel 7 with spin v and with a channel momentum ko-y determined 
by setting the channel energy equal to the system energy E: 



Ej{koj) 



E 



E^{p) = = Ei{p) + E2{p) = 



p^ 



p^ 



(16) 
(17) 



where the subscripts 1 and 2 refer to the specific particles in channel 7. 

The solution T of the Lippmann-Schwinger equation (^5|) automatically includes all iterations (ladder graphs) of 
the potentials but ignores crossed meson lines — which we assume to be small . The derivation of the potentials 
Vfja proceed via several step s. Fi rst in § [II A wc specify the coupled-channel Lippmann-Schwinger equation in the 
partial wave basis. Next in § 111 B we convert the Lagrangian (||) and (|9|) in to a Hamiltonian, and separa te off pieces 

j) and the contact potentials 



B we convert the Lagrangian (g|) and 
which produce interactions. Then we obtain the resonant potentials (§ 
the Foch-space matrix elements: 



HID) from 



13a 



'\HAa) 



Bo = S01,D03,P13 



]Hs\B^) 



1 



E~ M 



(BolHJa) 



Ba 



(18) 



The first term in dlq) describes direct scattering via an elementary quark transition in the contact interaction as 
illustrated in Fig. The Bq sum in ( p^ ) is over the processes, each second order in the Hg Hamiltonian illustrated 
in Fig. in which there are elementary resonance intermediate states as shown in Fig. ^ Resonances with bare 
masses Msoi, Mpis, and M^oa — but no widths — are thus built into the potential (as opposed to generated by it). 
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However, when the potential gets used as the driving term in the Lippmann-Schwinger equation (|5|), these elementary 
resonances get "dressed" , that is, the resulting T matrix contains resonances at shifted masses and with finite widths. 
Furthermore, we treat the elementary masses as adjustable parameters and let the data determine the best values. 
When the fitting is complete we then have dressed resonance energies and widths as determined from the behavior 
of the T matrix, as well as best-fit values for the bare masses. As we shall see, there is no guarantee that the bare 
masses and resonance energies are close. In addition, the potential derived from the contact interaction He generates 
"two-body" or "composite resonances" in T which are not explicit in the potential and thus not elementary. 



A. Partial Wave Lippmann-Schwinger Equation 

We assume the standard spin x i partial wave expansions for the T matrix 



T^«(k', /i'; k, = 2^ ^ [{I + l)Tl+{k', k) + IT'f^ik', k) 



X j dnk.dnkY{^,{k')Yi,^(k)TfU'^,tx'-'k,y) 

-e**'± sin(5;± 



1 ,.1 



(19) 



(20) 
(21) 



where x — cosd^k/, n — k x k', and T'* is shorthand for '^2. The partial wave expansion for the potential 
VaQ,(k', /i'; k, /i) has the same form as that for T^Q(k', /i'; k, /i) (p9|), so, for example, the projection for the potential 



1 ,,.1 1 ,.1 

im' ; ■ 

X / dnk'dnkY;„Ak')Yirn{k)Vf3^{k' , fi';k, ^i) 



(22) 



After substitution of the partial wave expansions of T and V into (|l|), we obtain the coupled one-dimension integral 
equations: 



9 roo 

TU^'^k)^Vl{k\k) + -Y, / dp 



pX',{k',p)T\=^{p,k) 



E + ie- E-y{k) 



(23) 



We indicated in Eq. ( [l4| ) that even for zero kinetic energy (the 1432 MeV threshold), a KN pair couples to nearby 
strangeness —1 charge channels. To use the derived potentials, we evaluate their matrix elements between charge 
states which in turn are expanded in the isospin states of those channels: 



K°n) = V2 

Y.-TT+) 
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0,0) + 




1,0) 








0,0) + 




1,0) 






V3 


0,0)' - 
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1,0)' H 


1 


2,0) 


-1 


0,0)' 


71 


H 




2,0) 


V3 


0,0)' + 


1,0)' H 


Ve 


2,0) 



(24) 



11,0)' 



Here the prime and double prime distinguish states with the same isospin but in different particle channels. Because 
our initial KN channel has no isospin-2 component, the |2,0) state does not couple to it and so we eliminate |2,0) 
and renormalize the rest. Using just two isospin channels for the three Stt states reduces computing times by 50%. 
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In order to include some isospin breaking, we do our calculations in the charge basis with physical masses for the 
particles and use these relations to transform between isospin and charge matrix elements. 

To solve the Lippmann-Schwinger equation ( p3| ) we have extended the Haftel-Tabakin technique (details in 
Appendix ^). First, rather than work with only the principal part we also keep the delta function (imaginary) part. 
Also, in handling coupled channels we reorganize the storage of matrices and instead of solving (|l^) we solve 

[G-^ -V][GT]^[V] (25) 

We now can use Gaussian elimination on the symmetric [G^^ — V] as opposed to the inversion needed to solve (|lj). 
Eliminating inversion saves ~ 30% in time and utilizing the symmetry of \V] and [G^^ — V] saves ^ 60%. 

The Coulomb force is included exactly in all our bound state (kaonic hydrogen) calculations, and of course it is 
crucial there. Even though we have the theoretical tools to include the Coulomb force in scattering, we do not include 
it in the total cross section calculations because the issue is too confused by different bubble-chamber experiments 
using different techniques for its removal during analysis. We do include some charge-symmetry-breaking effects by 
using physical masses, but otherwise we assume isospin is a good symmetry at the field theory and potential level 
and only break it by using physical masses in the energies for the kaon and lambda channels. 



B. Hamiltonian 



We deduce the Hamiltonian by following the canonical procedures starting with the energy-momentum tensor T^'^ 
2lPl: 

dC dC dC 

H = j (fx T°"{x) = Ho + Hs+ (27) 
Ho = Hmit + Hm (28) 
Hs = -J d'x ^qjt^j^XqdJ (29) 

He = Hot + Hcs 
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d'x jp<irX ■ (0 X do(b)q + I d'x -j^q > ^ • (0 X d,<ly)q (30) 



Here Hmit is the free bag Hamiltonian, Hm is the free meson Hamiltonian, Hs is the s-channel interaction (Fig. |l|A), 
and {Hct, Hcs) are the time and space derivative parts of the contact or four-point interaction (Fig. ^). We convert 
Hq into the Foch-space form Hq, by projecting onto the space of colorless baryons and expressing the meson field 
in terms of annihilation and creation operators wST 



Ho^ iBo\Ho\B'o) B'o (31) 

Bo ,B'g 

= E ^"^Bo + B^Bo + J2 [d^ku;kalik)a,{-k) (32) 

-Bo i 

Here Bq is the annihilation operator for a three-quark bag of type Bq, \Bq) is the bare baryon state, to_Bo is the MIT 
bare bag mass, and ujk = -^/ro^J^fA?^ is the energy of a free meson of momentum k. 



C. Resonance Potentials 



As discussed in §111 and illustrated by the Feynman diagram in Fig. g, for each meson-baryon channel a we identify 
one of the three resonance states Bq = A5oi(1405), Epi3(1385) or Ai303(1520) to include. As indicated in equation 
(Eq), we build potentials containing resonances within them by evaluating this Feynman diagram with the proper 



initial and final state quark wavefunctions. In Appendix ^ we simplify the Hamiltonian Hs which generates these 
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resonances, and in Appendices ^[^we evaluate the contribution to ( [T^ ) from each of the three resonances. After the 
partial wave projection there results the potentials 

{S01),ljfu i,\ - X X "Asoi/3(fe')'"Asoia(fc) ^0^^ 

32f^TTy/Uk'UJk E ~ Msol 
X {iB'is'^lM'hi' \IIb'Im' ; 00) (/sis; lMiM\lBlM;00) (33) 



X {lB'iB';lM'iM'\lB'lM';00){lBiB;lMiM\lBlM]00) (34) 
'x{lB'iB';lM'iM'\lB'lM'; W){lBiB;lMiM\lBlM; 10) (35) 

The vertex functions UBaa{k) for the 5 and resonances are seen in the appendices to be proportional to integrals over 
spherical Bessel functions (quark wavefunctions); for Spia that substitution is already made in (|3^). The potentials 
(P3|)-(|35|) are clearly separable and energy dependent. Further, each contains a pole at the real energy E = Mbo 
arising from the elementary resonances we have explicitly incorporated. As discussed in §111, for any value of Mbo 
the T matrices generated by using these potentials in the Lippmann-Schwinger equation will have poles at complex 
energies whose real parts differ from the Mb^'s, that is, the resonances get dressed and acquire widths. In addition, 
because the values for Mb^ are determined by the data fitting, their best fit values will differ from those we use as 
input. 



D. Contact Potentials v^""'^ 

As shown in Fig. |i|B, the contact interaction (|3^) directly produces BM — > B'M' scattering arising from an 
elementary quark transition Because B' and B belong to the baryon octet, the quarks in the initial and final 
states must be in Is states. We break the contact potential into the two pieces v^'^*'> and w^'^"^ which are the matrix 
elements of the time-derivative and space-derivative parts of He'. 

vfa^i^'^f^'^^'f^) = (/3,k',^'|i7ct,cs|a,k,/i) (36) 

For S*- waves the evaluation is straightforward: 

(k, k') = (5 wave) (37) 

i;l''^^(k',k) = >^YailB'iB';lM'iM'\lB'lM';IO){lBiB; lMiM\lBlM;IO) f," ^ ' (38) 

U3l\k',k) = N^iiOk+LOk') / drr^ [j^ito.r) + jUto^r)] jo{kr)jo{k'r) (39) 

Jo 

The coupling constants A^^ follow from the SU(6) wavefunction of the baryon octet |2^, and are given in Table ||. 

The increased spin and angular momentum coupling in P and D wave makes the vertex functions more complicated. 
The time derivative part is: 

z;(^^)(k',k) = ^] ^f{B'\ [ d^x6,q,,U,,,\,{uk+uJk.)fe'^'^-'^'y^q^s\BY^ (40) 



We evaluate it by substituting the Is quark wavefunction (Bl), substituting the partial- wave expansion of the plane 
wave, and integrating over solid angles to obtain: 



4t^(kV',k/i) ^ S,,,-^^f^^ Y.^;^^^)YUk')Nl I drr'[f,M + jfMjMkrMk'r) 

>^pi{lB'iB'; lM'iM'\lB' Im'; IO){lBiB; Im-imIIbIm; 10) (41) 



X 

/ 
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The partial wave matrix elements follow from the definition (p2|): 

v^;:^^''{k',k) = rdrr^ij^co^r) + jfM]MkrMk'r) 

^ IQnf-'y/ujk'UJk Jo 

X '^{lB'iB';lM'iM'\lB'lM'\Ii)>^Ya{^BiB; ImimIIbIm ; H) (42) 
/ 

The derivation of the space derivative part w^^^'(k',k) is more complicated than the time derivative part and we 
outline in Appendix ^ how: 

4a ''^(^''^) ^ /I ^2'/= = / drrjo{uJsr)ii{u)sr)ii{kr)ji{k'r) 

X ^ {Ib' iB' ; Im' Uf \Ib'Im' ; «) Aj^ (/sis ; iMiA-AlBhi ; (43) 

AH2V6/(^+i)(2/+i)(-i)^+'+i{l ! •} = {'^'t2l;^:!;| (44) 

The contact potentials (^Sf), (p2[), and ( ^ ) all are manifestly separable, energy (w) dependent, and contain vertex 
factors arising from the quark wavefunctions. They do not contain elementary resonances, although we will find them 
to be strong enough to generate composite resonances. 



IV. COMPARISON WITH EXPERIMENT 



We test the CBM by seeing how well we can reproduce five different groups of data after adjusting seven parameters: 
the bag radius i?, the couphng constants (/^"°, /^"^, /^"^), and the bare masses (Msoi, Mpi3, Mdos)- In comparison 
to potential models with adjustable coupling constants and ranges, the CBM's SU{i) x 5f/(3) symmetry greatly 
reduces the number of parameters — especially since the bare mass values mainly affect the local position of resonance 
peaks. 



A. Scattering and Reaction Cross Sections 



The first and largest data group is 300 measurements of two-body scattering and reactions cross sections for 
(70 < Kiab < 513) MeV/c = (1435 < E^m < 1567 MeV) ||,|^-||. These data, shown in Figs. |-|, are predominantly 
from bubble chambers and contain rather large statistical and systematic errors. The extension of the cloudy bag 
model to higher partial waves is important for understanding these data since even for a K~ momentum as low as 
150 MeV/c, P-wave contributions to K~p E^tt* have been reported |p2[, and at 390 MeV/c there is the A(1520) 



(45) 
(46) 



D-wave resonance (the peak evident in the figures near 400 MeV/c) [E6| 
The cross sections are related to the computed T-matrix elements (\ 




dapa 

dn 



k 



= 167r'*/x^^Q— ^ |T/3a(kV',k/i)|' 



<j0c.^ fd^^^ iQ^P^pf^ajY. \{i + 1) \Tl,l{k', k)\\i |r^;(fc', k) 

'' 1=0 L 



Here the bar indicates a sum over final and average over initial spin states, and the fi^s are relativistic reduced 
channel masses defined in Eq. (F3). For the channels we use the integrated cross sections are: 



a{K-p^K-p) = 16nfilJ2 [i^ + l^^lt (k' , k)\^ + I \tI^ {k' , k) 

1=0 

, , 2 

a{K-p -> ROji) = lQn^ll^l2- x ^ [(/ + 1) |T^^(fc', k)\^ +1 \T!^-{{k' , k)\ 



(47) 
(48) 
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<j{K-p ^ = 16n^l3^ll-J2 



1=0 



a + 
a + i) 



1 

1 



1 



V3 



/=0 



V2 



(49) 
(50) 



(51) 
(52) 



B. The A(1405) Resonance 

The second data group we examine is seven values of the E+tt^ mass spectrum determined from the K^p — > 
E+TT^TT+TT^ reaction at 4.2 GeV/c kaon lab momentum by Hemingway |Q. These data, shown in Figs. |^||, have 
a strong A(1405) resonance signal and are a major constraint on the model's parameters. Earlier A(1405) data 
also exist, but Hemingway's are cleaner since they were obtained through the three-step process K~p — > S+7r~, 
E+ A(1405)7r+, A(1405) S^tt". We calculate the Ett mass spectrum with a Watson model [P6|j34| which 
assumes an 5'-wave resonance dominates the final state interaction. The number of events per unit energy interval is 
accordingly: 



dN 
'dE 



cx /cctsx^Stt fx fc^3 |T3"3(fc,A:)| 



(53) 



We fix the normalization by setting the area under the theoretical mass spectrum equal to the total number of events. 
Because this model is so simple, we do not expect a detailed reproduction of the spectrum. 



C. The S(1385) Resonance 

The E(1385) is seen as a P-wave tt'^A resonance below the KN threshold. Aguilar-Benitez and Salicio have 
detected it in the mass spectrum from K~p Att^tt+tt" at 4.2 GeV/c. Our calculation of that Att spectrum is 
similar to the Ett calculation, only now with the Watson final state interaction model adapted to a P-wave resonance: 



dE k ^ k 



n^ikM (54) 



Because the resonance rests upon a large background, and because we have no model for such a background, we extract 
the resonance piece after fitting the entire spectrum as a polynomial background plus a Breit-Wigner resonance : 

(M, r) = (1384, 34.8)MeV, 004,2 = (-77554, 15303, -7402)MeV, h = 158.7MeV-'* (56) 
The resulting resonance is shown in Fig. as the peak at 1385 MeV. 
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D. The K p Threshold Branching Ratios 

The third group of data we examine is the branching ratios of K^p reaction rates at threshold: 

K~p^Y.^TT^ ^ charged K^p^n^A 

T " 1^ VI _ 7 Rc = — ^ Tj — , Rn = -TTz — r (57) 

A p ^ Zj+tt iv p — > all K p ^ neutral 

The measured ratios p^-|39| are shown in Fig. ^ The attraction of ratios is that they are not affected by uncertainties in 
target size and beam normalization, and so their uncertainties are smaller than those of cross sections. The difficulty is 
that calculations of rates into charge channels at zero energy are sensitive to Coulomb and charge-symmetry-breaking 
effects — not all of which are included in our model — and so we cannot be sure of the level of agreement to expect. 



E. Data Fitting 



The parameters of our model determined by various fits to these data groups are given in Table IV . The comparisons 
with the cross section and mass data are in Figs.^|], and the experimental and theoretical branching ratios are in 
Fig. ^. The data groups have large differences in the number of data points, the type of data, and the accuracy of the 
error estimates. Accordingly, it was not clear to us how properly to perform a best fit to all groups simultaneously. 



Consequently we have tried a variety of fitting procedures (indicated in Table III) in which different groups of data 
were systematically included and excluded, and with different choices for the x'' weights for each. 

In Fig. H we see fit 1. It employs only S'- waves and has its parameters adjusted to fit the low energy (p < 250 
MeV/c, E < 1470 MeV) two-body scattering and reaction cross sections. This fit provides a good base for comparison 
with the results of other groups and with our latter results in which higher partial waves are included. Fit 2 is also 
an S'-wave fit to low energy data, but extends fit 1 by including the Stt mass spectrum with its A5oi(1405) resonance. 
The fit labeled TRIUMF (A) uses the CBM parameter set A given by Veit et al. and uses the same criteria as 



our fit 2. Fit 2 and TRIUMF produce somewhat different parameters (row 1 ijs 3 in Table IV) due to differences in 
procedures and data. 

We see in Fig. ^ how much better an agreement fit 2 provides with the Ett mass spectrum than fit 1, and in 
Table |^ the large extent to which inclusion of the Ett spectrum affects the model's parameters. In examining the 
predictions for the branching ratios in the top part of Fig. H, we notice that fit 1 provides close agreement with i?„ 
and Rc, and fairly close agreement with 7. However, fit 2, which also included the Ett mass spectrum in its fitting, 
provides noticeably worse agreement with 7. For contrast, in the bottom part of Fig. ^ we present the branching 
ratios calculated with the Schnick-Landau potential model We sec that the updated parameters determined by 
Tanaka and Suzuki provide excellent agreement with the ratios. 

In Fig. H we notice that all three S'-wave fits fail to obtain agreement with the data at higher energies, and in 
particular, they fail to reproduce the Ado3(1520) resonance signal at 400 MeV/c. This clearly shows the need for 
higher partial waves and additional resonances. 

Our next two fits are shown in Fig. m Fit 3 uses 5, P, and D waves, examines data at higher momenta, and includes 
the three mass spectra for the A5oi(1405), Epi3(1385), and Ai303(1520) resonances. To determine the influence of 
the branching ratio data and to deemphasize the Ett spectrum, in fit 4 we do not fit to branching ratios and manually 
(the "m" in Table |l|) set the weight factors to {w scattering, wt.it, WK-n) — (6, 1,6). 

It is interesting to note in Fig. that fits 3 and 4, which include S — D waves, do not provide as good agreement 
with the K~p {K~p^K'^n,h^iT^) elastic cross sections as the S'-wave fits 1 and 2. What is clear from the figure, 
however, is that fits 3 and 4 provide nearly perfect agreement with the K^p Ett cross sections; apparently the large 
number of high energy, small error, Ett data points dominate these fits. This analysis is confirmed by confirming that 
the fit with the high energy data have a lower per degree of freedom. 

We again see in the Ett mass spectrum of Fig. ^ that fit 3 (which was adjusted to fit the branching ratios) does 
not produce a good Ett mass spectrum. Yet by releasing the branching ratio constraint we obtain fit 4 — which agrees 
with the spectra quite well. We also see that both fits do provide fairly good agreement with the ttA mass spectrum. 
This spectrum shows the sub _^A^-threshold Epi3(1385) resonance and essentially determines the value for Mpia. 

It can be argued that because so many assumptions go into applying the Watson final state interaction model to fit 
the various mass spectra, we should not expect good agreement with those spectra, and so should not use the spectra 
in data fitting. For these reasons, in fits n7 and cl we used all partial waves but fit only the two-body scattering and 
reaction data; n7 has weights determined by the experimental errors, cl has equal weights so the small Ett, error bars 
at high energy would have less influence. Again we note that without the Ett mass spectrum constraint, very good 
agreement with the branching ratios is obtained (fit n7 in Fig. 0) even though the ratio data were not included in the 
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fit. Yet if the Stt mass spectra must be fit, then fit 4 is our best fit even though we see in Fig. |^ that it does not fit 
the branching ratios well. 

Because the ttA mass spectrum determines the value for Afpia, and because this resonance does not significantly 
affect the cross section data above the KN threshold, not fitting this spectrum means the value for Afpia were not 
determined in fits n7 and cl. This explains the negative values they have in Table The other resonance masses 
do affect the cross section data above the KN threshold, and reasonable values were obtained for them. 



V. THE SCATTERING AMPLITUDES 



Now that the model parameters are determined, for each parameter set we examine the scattering amplitudes as a 
function of energy. We shall see that these amplitudes often have such a rapid energy dependence, especially near the 
KN threshold, that an examination of the full energy dependence is more significant that just the threshold values 
(scattering lengths). Values at that one energy are too sensitive to details such as channel mass values |^ . 

We start in Fig. ^ with the scattering amplitudes for KN — > KN in the Dq^ channel and for the ttA tt A in P13 . 
As we move down in energy we see vertical bars which signal the KN and Stt thresholds (the ttA threshold at 1252 
MeV is not visible). We also see cusps in / at channel openings and that Im / — but not Re/ — vanishes when all 
channels are closed. The KN resonance A£)03(1520) and the ttA resonance Epi3(1385) are clearly signaled by Re/ 
changing sign and Im / peaking at approximately the energy corresponding to the resonance mass. The bare versions 



of these resonances were built into our effective potentials via (18), yet this does not guarantee that they will remain 
after their dressing by the Lippmann-Schwinger equation and the adjustment of parameters in the data fitting. 

In Fig. H we give the / = 0, iS-wave KN scattering amplitude for fits 1-4 and TRIUMF-A (since fit 2 is an updated 
TRIUMF-A, the two are very similar). We see that fit I's Re/ changes sign once below the KN threshold (right-most 
arrow in Fig. ^ and once again slightly above the Ett threshold (left-most arrow). All the fits' Re/'s pass through 
or very close to zero above and below the KN threshold, but only fit I's Re/' has an extra sign change far below. 
This is fascinating, for the ability of the model to reproduce the experimental IS* strong interaction shift of kaonic 
hydrogen depends on the magnitude and sign of Re/ at threshold, and clearly there are several sign changes here. 

The imaginary parts of / (lower part of figure) for fits 1-3 show a shoulder right below the KN threshold at ~ 1440 
MeV, while fit 4 shows a peak at ^ 1425 MeV. None of these behaviors for Re/ and Im / are the classical resonance 
signal expected at 1405 MeV for the A5oi(1405). Although there appears to be a resonance signal somehow mixed 
into the KN threshold behavior. Fits 1 and 3 do show what appears to be distinctive resonance signals in both Re/ 



and Im /, but they appear way down in energy near the Ett threshold of 1325 MeV! Clearly, since Table IV shows the 
fitted Msoi > 1550 MeV, the cloudy bag model does not give the Asoi(1405) as a simple bag state of mass ~ 1405 
MeV. 

In order to unravel some of the mysterious resonance behavior of the ^oi KN scattering amplitude, we have solved 
for the complex energies at which the scattering amplitudes have poles. If the pole is on the second (unphysical) 
energy sheet near the real positive axis, it will produce an observable, narrow resonance along the positive real energy 
axis with a width related to the imaginary part of the complex energy by 

r ~ -2Im E (58) 

Other singularities or channel openings can interfere with this picture. If the pole is far away from the real energy 
axis or some other singularity gets in the way, its experimental and dynamical significance is less. In Appendix ^ we 
show that the formal solution of the Lippmann-Schwinger equation implies that the pole positions are solutions of 
the equation 

det(l - VeGe) = (59) 



We find numerically the complex energies which solve (59) 



In Fig. y we have plots of the imaginary parts of the KN S wave scattering amplitudes of fits 2 (upper part of 
figure) and 3 (lower part of figure) as functions of complex energy. As found in the explicit solution of (|59|), and 
evident in the figure, there are two poles in the 5*01 channel at the complex energy listed in Table ^ For fit 2 we 
note that the high energy pole is actually 17 MeV above the KN threshold (the tick mark to the right of 1400 MeV 
in Fig. and 31 MeV along the negative imaginary axis. Fit 2's low energy pole is seen to be 80 MeV below the 
KN threshold and far from the real energy axis. In contrast, fit 1 and fit 3 (in the lower part of the figure) have low 
energy poles very close to the real energy axis. They are the causes of the striking energy dependences near the Stt 
threshold seen in Fig. ||. 

In contrast to these CBM fits, we note in Table ^ that the potential models of Alberg et al. and Schnick and 
Landau [p9| have only one S-wave pole, and it is closer to the tabulated A(1405) energy. Apparently in the potential 
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model, the A(1405) is a composite resonance with a single pole close to the tabulated energy, while in the cloudy bag 
model, the A(1405) is not elementary and not simple, a conclusion drawn previously from a number of viewpoints 



|2|,^J34[|4C|,^. 

The T matrices' pole positions for the P and D wave amplitudes are also listed in Table 0. These positions 
correspond to a I]pi3(1395) dressed mass of 1385 MeV and width of 24 MeV, and a A£)03(1520) dressed mass of 
1514 MeV and width of 12 MeV. The values for the best-fit bare masses are Mpis ~ 1419MeV and Mdos - 1552 
MeV. This means that in both cases the renormalization by the contact interaction and higher order scatterings shift 
the masses downwards in energy by 31 MeV!. The tabulated masses and widths of (1383.7 ± 1.0, 36 ± 5) MeV 
and (1519.5 ± 1.0, 15.6 ± 1.0) MeV are in excellent agreement with the pole positions — especially since the two are 
expected to differ somewhat when the pole is not close to the real energy axis. The tabulated half widths of the 
5^pi3(1395) and A£)03(1520) are not in as close agreement with the imaginary parts of the pole energies (18±3 vs 12, 
8 ± 1.0 lis 6), but this is expected since the widths arise completely from renormalization, and this is a broad and 
nonsymmetric resonance. 



A. Kaonic Hydrogen 



There are three measurements of the width F and strong interaction shift relative to Bohr energy, e = —{E — Eb), of 
the IS* level in kaonic hydrogen |l^]- jl^. Although these measurements should be a good test of the KN interaction 
slightly below threshold, the uncertainties in the measurements make conclusions difficult (we see the statistical 
uncertainty in Fig. |l^, but not the apparent systematic uncertainties). Since at present there is distrust of a theory if 
it agrees with these experiments, we eagerly await the new experiment in progress with its promise of high precision 
and high accuracy Q. 

In our calculation of the kaonic hydrogen state, we identify the complex bound state energy as the T matrix pole 
energy for the combined Coulomb-plus-nuclear force problem , and solve ( |59| ) for them. Our predictions are shown 
in Fig. For all the CBM fits the calculated width F is acceptable or slightly large, yet the shift e is of opposite 
sign to the data (all experimental shifts are to the more bound). If the data are correct, then this implies that the 
the real part of the K^p scattering amplitude (a sum of / = and 1 = 1 amplitudes) is positive at K~p threshold 
even though the / = piece dips below zero at threshold. The potential model of Schnick and Landau |jl^ or its 
update by Tanaka and Suzuki EG] have this property, that is, they agree in sign with the data. Although the sign 
changes in some of the cloudy bag model's Re/501 are similar to those of the potential models, none of our CBM fits 
have the proper combination of / = and 7=1 strengths to keep Kef{K~p) > 0. Although we have not done it, we 
suspect that if this condition were made a requirement of the search, such a solution would be found (presumably at 
the expense of other data). 



VI. SUMMARY AND CONCLUSIONS 



We have extended the SU (3) cloudy bag model for the coupled {KN, Stt, Att) system from S to _D-waves and thus 
to much higher energies. The model has elementary quarks inside a bag as well as an SU (3) meson field inside and 
outside the bag. While not derived from quarks, the meson field restore chiral symmetry. The model Hamiltonian 
contains a contact interaction which generates direct meson-baryon scattering as well as an s-channel interaction used 
to include elementary A5oi(1405), Spi3(1520), and Ai303(1385) resonances. We derived effective, energy-dependent 
separable potentials for use in the Lippmann-Schwinger equation as well as some new numerical approaches to solve 
the coupled Lippmann-Schwinger equation. The parameters of the model were determined after extensive fits to 
various scattering, reaction, branching ratio, and mass spectra data. 

Our fitting and subsequent analyses indicate that the 2(1395) and A(1520) are well described as elementary reso- 
nances. Although their pole positions are renormalized by the contact interaction and higher order scatterings from 
the Lippmann-Schwinger equation to 31 MeV below their bare masses, the resonance energies (as determined by the 
complex energy poles of the T matrix) are within 2 and 4 MeV respectively of the tabulated values. 

The cloudy bag model's description of the A(1405) S-wa.ve resonance is less simple. While we agree with the 
previous conclusion that the state is not an elementary, three-quark s-channel resonance ||l^,|l^,^,^,Q , we have 
also found that it is less simple than the quasi-bound KN, single pole state produced by potential models ||l^,Q. In 
particular, there are two poles present in this channel, with the resonant behavior near threshold arising from a pole 
above threshold interfering with the threshold cusp. 

Although the two sign changes in some of the cloudy bag model's KN scattering amplitude near threshold are quite 
similar to those of potential models which agree in sign with the strong interaction shift in kaonic hydrogen, none of 
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our CBM fits have quite the proper combination of / = and 1 = 1 strengths to keep Iief{K~p) > 0. Agreement 
with the shift could be required as part of the fitting procedure, but we suspect that this is best left for a time when 
more acceptable data are available [Q. 

In a general sense we conclude that the CBM with S, P, and D waves is able to reproduce the K^p scattering 
and reaction cross sections from 70 — > 513 MeV/c (1435 1567 MeV) and either the K^p — > (Stttttt, Atttttt) mass 
spectra or the branching ratios 

K^p^Ti^TT^ ^ ^ charged ^ K^p^n^A . ^, 

^ = iF^ — ^ 11 — ' = iF^ r~r f-^^) 

K p ^ L^TT iv p — !■ all K p ^ neutral 

Fitting all three together appear too much to ask from such a simple model. 

We note that the fitted parameters appear reasonable. The average bag radius R ~ 1.22 ± 0.14 fm is about 0.1 fm 
larger than the value found by Veit et al. jlj], but within the range 1.5 > i? > 1.0 fm given by Guidry (the range 
is for simple to refined models, and our model with massless quarks is simple). If we restrict our fit to low energy 
scattering and mass spectra data we obtain R = 0.95fm, which is not big and indicates that larger R values arise from 
the effort to fit the high energy data. The fitted meson decay constants (/■^"°, /^^^j — (99, 88, 91) MeV appear 

quite close to accepted Q values {fx, fw) — (112, 93) MeV, especially since SU{3) predicts = fx- 

In conclusion, we believe the cloudy bag model has success in reproducing much data — but clearly with limits. The 
difficulties may arise from trying to reproduce complicated energy dependences with too simple a model. We have kept 
only linear terms in the meson field and have assumed a rigid and spherical bag of one size for all baryons. The bag 
does not recoil and the point-like mesons do not arise from quarks. We expect this to cause difficulties at the higher 
energies and higher momentum transfers. We have used the MIT bag wavefunctions for massless quarks in a square 
well, and so have form factors proportional to spherical Bessel functions. This too is rather restrictive. We have looked 
at zero kinetic energy branching ratios, but have included isospin-breaking effects only at the Lippmann-Schwinger 
equation level not at the Hamiltonian level. Clearly, all these effects may be important in a more sophisticated study of 
branching ratios and of the kaonic hydrogen level shifts. Further improvements appear worthwhile as do applications 
of the model within a nuclear environment. 
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APPENDIX A: REDUCTION OF Hs 

We convert Hg (v% into a more calculable form by partial integration: 



H. = - d 



(Al) 



Application of the Dirac equation for the bag (H) shows the last term to vanish. The space derivative in the first term is 
converted to a surface integration over a infinitely large surface, which in turn vanishes since the quark field is confined 
within the bag. The linear boundary condition (^ and the surface delta function imply —{dn9y)q^^^^\q ■ (j)/2f = 
— ■^Sgq^y^Xq ■ (j) and this yields: 



Hs 



g75A • q<j)Ss - ^ 9o(g7"75A ■ 



(A2) 



We use this Hg in the derivation of the potential to follow where we separate off a vertex function V in the Foch-space 
Hamiltonian for BAI ^ B'r,: 



Vo^{k) = B'„UBo\Hs\a) Bo 



Hs = Y.j [^o.(k)a,(k) + yj,(k)aj(k)] 



(A3) 
(A4) 
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APPENDIX B: BM ^ Asoi(1405) POTENTIAL 



The quark model configuration for the bare A5oi(1405) is one u, d, and s quark in an SU(3) flavor singlet, with one 
lpi/2 quark and two Is quarks. The transition BM — > A501 thus has one quark absorbing a meson and being excited 
from the Is to lpi/2 state. The quark wave functions are: 



qfAr.t) = 4^ ( .A'^f , ) e— =*0(i?~r) (Bl) 



Here x is the spin-flavor wave function of the quark, jo and ji are spherical Bessel functions, (uJsjUJpi) ~ 
(2.04/i?, 3.81/i?) are the energies of Is and lpi/2 states, and the iV's are normalization constants Efl]. The 
BM ^ A5oi(1405) vertex function needed in the Foch-space Hamiltonian Hg ( |A4| ) and the form factor for the 
Hamiltonian are accordingly: 

Vi'°'\k) ^ A\/^v^,,,^ik) Bo (B4) 
= {^SQi\Hs\a) (B5) 



NsNpi { 2R^jo (w.i?)jo i^piR)jo (kR) 



2/ 7(2^0^ 

+ (wfc +uJs- LUpi) J drr"^ [jo{i^sr)jo{i^pir) + ji{uOsr)ji{uJpir)]jo{kr)^ (B6) 

The spin-flavor matrix element '^^{Ksoi\K\B)'^^ is evaluated with the Wigner-Eckart theorem: 

^■^(Asoi|A,|B)^/ = A^(/B^s;/Af^A^|/B/M;00) (B7) 

where the bracket on the right-hand side is a Clebsch-Gordan coefficient. The coupling constants are calculated 
using the SU(6) quark wave function [|l2|.^ and are listed in Table ^ The potential u^'^'^-'^' (k', /i'; k, /i) corresponds 
to meson-baryon scattering through an intermediate A^oi state, FigTM: 



V 



(SOI) 

_ 5titi'^p^t{lB'iB';lM''i'M'\lB'lM']^^) {lBiB\lMiM\lBlM\ 00) UAsai0{k)uAsoia{k) 

E-Msol 64/27r3 V^^^i^ 



\H,\Asoi)-e; ^^{Asvi\Hs\a) (B8) 

(B9) 



The partial-wave matrix projection of this potential via d23) yields (35) 



APPENDIX C: BM ^ Ai303(1520) POTENTIAL 

The A£)03(15 20) has one u, d, and s quark in an SU{3) singlet with total spin 3/2 and wavefunction Is^ lp3/2- The 
KN —^ A^ios transition thus has one quark absorbing the antikaon and changing from Is to the IP3/2 (flavor may 



also change). The vertex function for the Foch-space Hamiltonian (A4) and the IP3/2 quark wave function are: 



y(^°^'(k) = At , / ''(^^0^\^sqpl,Kqs + {COS ^p,)qpl'l,K<ls\BYf 

(r, t) - Np, ( . ) eiR - .)e--3*3;f| (6, 0) (C2) 

Here 3^/^^ is the spin-angle function, A'p3 is the IP3/2 normalization constant and ajp^ « 3.20/i? is the energy 
level of the 1^3/2 state. We express the spin-angle function in terms of spherical harmonics and Pauli spinors, 
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substitute \Ado3}^^ = X^mil"^' ~ "^I^^S i-^)^inil^Soi)''''' Ado3 spin-flavor wavefunction (the same as 

Asoi) to obtain: 

^^(Asoi(m " m)\X,a^\Biiy)yf = -Vsxti^:^; lq\^l; ^fx ~ m){lBiB; hiiMllBhi; 00) (C3) 

1 1 47rNsNp3 v^v^ 1 ,13,^,-, 

^A.c3a(M, ^) = -^^-j==== l^ ^ (-l)''(lm; - m|l-; 

^ V V ) rnq LAI 

X 'f{Asol{^i-m)\X,ag\B{,^)yf ( (fxYi-q{r)Y*^{f)YLM{r)jL{kr) 



X {(5s[ji(ws7-)ji(wp3r) jo(uJsr)j2{uJp3r)] 

+ 6'„(tJs + c^fc - c^p3)[ji(ws?')ji(wp3?') - jQ{.'^sr)i2{<~^pZ'r)]} (C4) 

= —l=^^(-ir*^V2rTTy£M(fc)4';lo(fc)(^o; 101^1; 

•' ^ ^ LM mq 

X ^^(Asoi(m - m)\\,a^\B{v)Yl {Im- i/i - m|li; ^/i)(LM; l-m|il; Iq) (C5) 
'^So3o(^) = -NsNp^ [2R^j^{LOsR)3i[LOp3R)iL{kR) 

+ {uJs + Wfe - Wps) / drr^jL{kr)[ji{usr)ji{ujpzr) - jo(ws?')j2(wp3?')] (C6) 

JO 

The coupUng constants A„ are the same as in the BM A501 case (Table and after lengthy algebraic manipulation 
of the Clebsch-Gordan coefficients and 6j symbols, we obtain the vertex function: 

?^AD03a(k;A«,i^) = , J. °, — {Ibib ; hiihAlBlM ] 00) 



X 5^5^(-l)"z^V2mr;,,(fc)(L0;10|Ll;10)^i^„3Jfc) 

X (Im; - TO|li; ^n){LM; l-m\Ll; Iq) (^i/; lizl^l; - m) 
= {IB^B■, Imim\IbIm; 00) (1^; 2^^ - H^2; \^.)Y;(,_,){k) (C7) 

The separable potential corresponds to meson-baryon scattering through the intermediate h-om, Fig- |[ i-G- 

using Auos for Bq in (^: 

^;^^°^)(k^k) = 5](/3(M')l^.|AD03(^));^3^(Az,o3(^)|i^.|«(M)) 

_ A^A^MAD03/?(fcO'"AD03a(fe) {lB'iB';lM'iM'\lB'lM';00) {IbIb; iMiMjlBlM^OO) 



V 

The partial wave projection of Tjf^"^) follows pages of algebra and yields (0). 

APPENDIX D: BM ^ Epi3(1385) POTENTIAL 

The bare Epi3(1385) is a single u, d, and s quarks in a 10 representation of SU(3), with all quarks in the Is state. 
The transition BM — > Spis thus has one quark changing its flavor-spin state after absorbing a meson. The vertex 



function and form factor for the Foch-space Hamiltonian (A4) are accordingly: 
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K,f>(k)=E°;t^Sp,3a(k)So (Dl) 

i;sp,3„(k) = (I]pi3|ir.|a)So (D2) 
1 /" d'^r 

= — / ^===^=e*'^'"(I]pi3|i(5sgis75Aigis + i6'«Wfe9is7°75Aigis|S) (D3) 

= 2;i^(-2) y ^====e^''-'-JoKr)ji(^,r)5//(Spi3|A,:^ • f\BY^ (D4) 



1 A^? 



X. TAT 



2/ 4^ 7(2^0^ 



1 

t2 d2 



LM 

I An, 



TTT^ E n;(^)joKi?bi(c..i?)ji(fci?)^-^(Spi3|A,a,|i3)^/ (D5) 



To obtain the vertex function, we substitute ( Bl ) for the Is quark wave function, substitute the partial wave expansion 



of the plane wave, substitute the explicit form for cr • f, and use the Wigner-Eckart theorem: 

^^(Epi3(M)|A,a,|SH)^^ =A^(ij.;l(Z|;il^A^)(/BiB;/M*M|/B/M;10) (D7) 

—i)spN^li^ 1 13 

?^Spi3Q(k;Ai,i^) = o ^ (o'^' l^'lo^! 7;l^){lBiB-:lMiM\lBlM] 10) 

X nV-.) (fc)io(^.i?)jiKi?)ii(fci?) (D8) 
The calculated coupling constants A„ are given in Table The potential w^-^i^) (j^/^ ^) 

is similar to ^ 

now with an intermediate Spi3 (Fig. 



4™)(k',k) = ^(/3(A.')|gs|Spi3(^)) p ^ . (Spi3(^)|g.|aM) (D9) 



127r2/2^/ZJ^ S-A/P13 

^ ^ ~1 131 13 

^ H^iV-M) W^i('^-M')(^')(2A*; - rnu\-l; i^^) ii^l^-' - mu'\-l; -v) 

X {lB'iB';lM'iM'\lB'lM'-AO){lBiB;lMiM\lBlM; 10) (DIO) 

Its partial wave projection via ( p2[ ) yields (^5|). 

APPENDIX E: CONTACT POTENTIALS 

The space-derivative part of the contact potential is: 

= / drdr!,rjo(a;..r)ji(L^,r)^/(B'|/»'»jA,'T|S)^/ 
47r 7o 

e~*'^'-'"(r X k)e^'' '' + (r x k')e-''^' ""e*"' 



4(2^)3/2 VZJ^ 



(E2) 
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where we have substituted for the quark wave functions qis and qi^. Since r x k and r x k' are the angular momentum 
operators acting on e^^ '^ and e*'' '", substituting the partial wave expansions of the plane waves yields: 

<^)(k',k) = ^^'^1^ yLimm'Ck') rdrjoMnMMkrMk'r) (E3) 
A = 'f{B'iti')\f,,,,X,a\Biti)yf J dQrY;^,LYir,, 

= E(-l)' / dnrY:^,L,Yir,^f{B'ifj,')\U,,X,a.,\Bit,)yf 

f{B'\\u,,x,cT\\B^f + mi + 1) { I f ]] 



X (^M; lm\^l;jm' + lm'\^l;jm' + fi') (E4) 

where we have used the tensor notation for a-q and Lq as in ( |D(]| ), the Wigner-Eckart theorem, and many manipu- 
lations. The calculated coupling constants A^'^ are given in Table |l[ After substituting all these relations back, we 
obtain: 

t;£^(k',k) = -N^ J2 ^"^iilBiB-, Imim\IbIm; IO){lBnB'; hriM'llB'hp; 10) 

I .jM.lmra' 

,1 , ,1, ,1 , , „i, 



X (-/x; lm| -l■]M){-^i'■ Im' j M)YC„,{k)Yi^, {k')^Ql{l + 1){21 + I) 
y,(-\\3+i+h f"^ drrjo{LO,r)n{^,r)ji{kr)ji{k'r) \ \ \ ll 

^ ^ io I ; ? J- / ^^^^ 

The equations for w*^'^''-' and w*^"^*-* have the same form as equations (2.8) and (2.9) in Ref. but the coupling 

constants (spin-flavor matrix elements) in Table || are different. Finding the partial wave matrix elements of the 
potential v^'^^^ is more complicated. We start with the definition ([l9|), apply the orthogonal relations of the spherical 
harmonics, the properties of the Clebsch-Gordan coefficients, and lengthy manipulations to obtain (Esl). 



APPENDIX F: EXTENDED HAFTEL-TABAKIN TECHNIQUE 



We give here our extension of the Haftel-Tabakin technique |2J] for solving the Lippmann-Schwinger equation (|2^). 
For each channel 7 the integrand in ( ^ ) contains an integrable singularity at the "on-shell" channel momentum 
/co7 = ^0 defined in ([l7|). To permit numerical integration this singularity is removed by subtracting a term from the 
integrand which leaves the integrand nonsingular ||24|, and then adding in the integral of the subtracted term: 



I' 

Jo 



dp 



p'V{k',p)T{p,k) 
E - E{p) -f ie 



dp 



p-'V{k',p)T{p, k) 2fiklV{k', ko)T{ko, k) 
E - E{p) + ie 



fr2 

"-0 



2^jikf,V{k',ko)T{ka,k) , -.,^2_p2 

p'V{k',p)T{p, k) 7rMo^(fc', fco)T(fco, k) 



dp 



E - E{p) 



k^ 
ftp 



-TrifikQV{k\ A:o)r(fco, k) 



(Fl) 



(F2) 



where the RHS of Eq. (Fl) is evaluated analytically and where /i = /i^ is the relativistic "reduced mass" for channel 



7: 



M7 = o 



1 dp^ 



2 dE^ 



E-y E 



Ei{ko^)Ei{koy) 
Ei{ko^) + Ei{koj) 



(F3) 



We solve the integral equations (|2^) by approximating the integrals as sums over N Gaussian grid points {pi\i — 1, N} 
with weights {wi\i = 1,N}: 
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Tpaik', k) = Vi3a{k', fc) - X! '^'ifJ"iko^yi3f{k', koj)Tja{koj, k) 

7 

pfVp.y{k',p^)T^a{p„ k) 2fi^k^^Vi3j{k' , koj)T^a{ko.y, k) 



TV 



E-E^{p,) 



We convert (F4) to the set of linear equations, 



T, 



Vr 



nj3.ma 



EE 

7—1 i—1 



Vr 



np.ij 



k^ 



(F4) 



(F5) 



by defining the supematrix elements for G, V , and T 

Vf3a{Pn,Prn) 
5 Pm } 
V(3a{Pn, koa) 
Vf3a{kQf3, koa) 

_ \ 2WnPll{TT{E ^ 



-2^*7 E 
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(n 
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in 

E^iPn)) 



IN, 
N + 
IN, 
N + 



m 
m 
m 
m ■ 



IN) 
IN) 
N+1) 
N+1) 



{^Wikl^ 



/iAk'o, 



pf)) ~ ikQ^2fi^ 



{n = lN) 
{n = N + 1) 



(F6) 



(F7) 



Here Tnp.ma are the N^{N + 1)^ unknowns for Nc channels and N Gaussian grid points (the +1 arising from the 
on-shell point). The potential V is taken as the sum of the 5 terms in (|l^), Nc = 5 for the five independent isospin 
channels in (E4), and we use 24 or 32 Gaussian grid points in each channel. We then solve the N^iN + 1)^ coupled 
equations (|F5|) by rearranging the matrix equation: 



[T] = [V] + [VG][T] 
=> [l-VG][T]^[V 
^ [T]^[1-VG]'^^[V] 



(F8) 
(F9) 
(FIO) 



Because the solution for T is numerically intensive, we have made a number of improv emen ts to our former technique. 
First, while we previously determined the actual inverse [1 — T^G]"^ and evaluated ( FlC ), we now solve (^) for T 
using Gaussian elimination. This is ~ 30% faster than matrix inversion. Next we obtain an even greater savings by 
utilizing the symmetry of the potential V and the diagonal nature of the Green's function G to reformulate ( |F9| ) as 

[G-^ - V][GT] ^ [V] (Fll) 

Since [G^^ — V] is symmetric, we solve a symmetric linear system of equations for [GT] — which is much faster than 
a general system. From this [GT] we easily obtain the value of [T] since [G] is a diagonal and its inverse is quick to 
compute. 

If there are some elements of G which vanish, the reformulation is a little more complicated. For example, we 
assume the last element of G vanishes while none of the others do (if it is not the last element, we can always 
rearrange columns and rows to make it last). The linear system is now: 

(F12) 

Here v is an {M — 1) x {M — 1) symmetric matrix, g is a {M — 1) x {M — 1) diagonal matrix, t is a {M — 1) x {M — 1) 
matrix, Vr and tr are (M — 1) x 1 row matrices, Vc and tc ar e (M — 1) x 1 column matrices, and vm and tM arc the 
lower right-hand corner elements of V and T. We split ([F12|) into the equations: 



l-vg 


0' 




' t tc' 




V Vc 


-Vrg 


1 




tr tM 




Vr Vm 



(1 - vg)t = V 

(1 - vg)tc = Vc 

— Vrgt + tr = Vr 
-Vrgtc + tM = Vm 



(F13) 
(F14) 
(F15) 
(F16) 



and solve for t and tc with the symmetric matrix v (this is fast). We then solve for tr and tM if needed. If there 
are several elements in G which vanish, we follow the same procedure recursively until we are left with a symmetric 
system to solve. In practice, the direct LU decomposition method took 1214 second on an IBM RS/6000 Model-530, 
while the method using symmetry took 493 seconds. 
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FIG. 1. Fcynmau diagrams for meson interactions in the cloudy bag model, (A): via the s-channel Hamiltonian Hs, (B): via 
the contact Hamiltonian He- Dashed lines represent mesons and solid lines fermions (quarks or baryons). 
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FIG. 2. The Feynman diagram used to generate potentials corresponding to intermediate resonance excitation via the 
Hamiltonian Hs- 



FIG. 3. The K~p scattering cross sections, reaction cross sections, and Ett mass spectrum calculated with the S wave fits 
1, 2, and Triumf (A). As is evident, fit 1 was not adjusted to fit the mass spectrum. The cusps right below 100 MeV/c arise 
from the opening of the K*^ N channel, and the peaks near 400 MeV/c arise from the A(1520) _D-wave resonance (not included 
in these 5'-wave fits). The cross section data are o: Cibrowski et al. □: Kim j28|, o: Mast et al. js^. A: Sakitt et al. [ p7[ , 
<: Watson et al. [^, V- Evans et al. >: Kittel et al. |2^, and x: Bangerter et al. js^. The Ett mass spectrum is from 
Hemingway p^ . 

FIG. 4. The K~p scattering cross sections, reaction cross sections, Ett mass spectrum, and tt^A" mass spectrum calculated 
with parameters of fits 3 and 4. Both fits include S, P, and D waves, but fit 4 was not adjusted to branching ratios (shown 
in Fig. The cross section data are the same as in Fig. ^ with the addition of the tt^A" mass spectrum extracted from the 
results of Aguilar-Benitez and Salicio [ ^ . 

FIG. 5. The K~p scattering cross sections, reaction cross sections and Ett mass spectrum calculated with parameters of fits 
n7 and cl. Both fits include 5", P, and D waves, but were not adjusted to mass spectra data. Fit cl used equal weights for the 
data in the minimum search. The cross sections data are same as in Fig. ^. 

FIG. 6. The threshold branching ratios for K~p to neutral states R„, to charged states Rc, and to charged Ett states 7. The 
comparisons in the top part of the figure indicate that CBM fits adjusted to the Ett mass spectrum tend not to agree with 
7. The comparisons in the bottom indicate that the original Schnick-Landau potential parameters did not provide good 
agreement for 7, but the update (solid curve) by Tanaka and Suzuki does. The data are from Humphrey et al. Tovee 
et al. Nowak et al. Q, and Goossens et al. psj] . 

FIG. 7. The P13, ttA and D03, KN scattering amplitudes of fits 3 and 4. The Im / in the lower halts of the figure peaks at 
the resonance energy for which the Re/ passes through zero in the upper part. The arrows show the Ett and KN thresholds. 

FIG. 8. The .^A'^ Soi scattering amplitudes for fits 1-4. The arrows show the Ett and KN thresholds. The resonance and 
threshold behaviors appear mixed. 

FIG. 9. The imaginary parts of the KN S-wave scattering amplitudes of fit 2 (upper part) and fit 3 (lower part) as a function 
of complex energy. The projection onto the complex energy plane shows the contours of the T matrix. The tick mark to the 
right of 1400 on the real energy axis is the KN threshold energy. 

FIG. 10. The shift and width due to the strong interaction of the IS level in kaonic hydrogen. The data are from Davies et 
al. |l^], Izycki et al. ||l^, and Bird et al. jl8|. The triangles are prediction of the different CBM fits. All experiments indicate 
a positive e while all CBM calculations indicate a negative shift (to the less bound). 
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TABLE I. Channel assignment for charge and isospin bases. 
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TABLE II. The coupling constants A^^ and A^ 
and p. 


for channel a, 
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TABLE III. Characteristics of different fits. 
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1.29 
0.95 
1.00 
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1.34 
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1623 


1421 
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TABLE IV. Parameters of different fits; R is in fm., others in MeV. 




Fit 










E{S) 






E{P) 






EiD) 


Vahw 


1 
2 
Tri 

|15I 

m 

3 

4 








(2, -27), (-104, -6- 10"") 
(17, -31), (-80, -76) 
(13, -12), (-95, -76) 
(-52, -39) 
(-14, -25) 
(5, -21), (-129, -4) 
(1.2, -13), (-114, -74) 






(-48, -12) 
(-46, -12) 






(81, -6) 
(83, -6) 



TABLE V. The pole positions in complex energy (in MeV) for S — D waves relative to the KN threshold at 1432 MeV. 
Numerical uncertainty is ~ ±2 MeV 
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